Quantifying turbulence induced segregation of inertial particles 
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Particles with density different from that of the advecting fluid cluster due to the different response 
of light/heavy particles to turbulent fluctuations. This study focuses on the quantitative charac- 
terization of the segregation of dilute poly-disperse inertial particles evolving in turbulent flow, as 
obtained from Direct Numerical Simulation of homogeneous isotropic turbulence. We introduce an 
indicator of segregation amongst particles of different inertia and/or size, from which a length scale 
r seg , quantifying the segregation degree between two particle types, is deduced. 

PACS numbers: 47.27.-i, 47.10.-g 



The ability of efficiently mixing transported substances 
is one of the most distinctive properties of turbulence, 
which is ubiquitous in geophysical and astrophysical flu- 
ids. New features appear when turbulent flows are 
seeded with finite-size particulate matter having den- 
sity p p different from the carrier fluid density pf. Due 
to inertia, measured by the Stokes time T = a 2 /(3/3is) 
(a being the particle radius and v the fluid viscosity; 
/3 = 3pf/(2pp + pf)), such particles detach from fluid 
parcels' paths and distribute inhomogencously [1, 2, 3]. 
Although this phenomenon of preferential concentration 
[4] has been known for a long time [1, 2], it continues 
to attract much attention (see [3, 5, 6, 7, 8, 9, 10, 11] 
and ref. therein). It is important for drag reduction by 
microbubbles [12], for the effects of microbubbles on the 
small scales of turbulence [13] , for aerosol physics which is 
critical for climatological models [14], or to understand 
the patchiness of chemical and biological agents in the 
oceans [15]. The key issue is the tendency of inertial par- 
ticles to form clusters with the consequent enhancement 
of the particle interaction rate. 

When having different particle types in the same flow 
(polydispersity), the respective particles probe different 
flow structures: light particles (f3 > 1, e.g., air bubbles 
in water) preferentially concentrate in high vorticity re- 
gions, while heavier ones (/?< 1, e.g. sand grains in water) 
are expelled by rotating regions. This leads to a segre- 
gation of the different particle types, which intuitively is 
characterized by some segregation length scale. An exam- 
ple of particle segregation is shown in Fig. 1, where snap- 
shots of light and heavy particles' positions are depicted. 
The segregation length depends on both the respective 
particle densities and Stokes numbers St = r /t^, which 
measure the particle response time r in units of the Kol- 
mogorov time T n (characterizing the smallest active time 
scale of turbulence). This Letter aims to systematically 
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FIG. 1: Slice 4OO77X4OO77X lOry of heavy = (red) and 
light j3 = 3 (blue) particle positions. From left to right 
St= 0.1, 1, 4.1. Data refer to the simulation at Re^ = 180. 

quantify the segregation length as a function of both the 
relative density (/3) and the Stokes number, which to- 
gether characterize the particle classes. 

To demonstrate our method, we consider a model of 
passively advected dilute (to neglect collisions) suspen- 
sions of particles in homogeneous, isotropic turbulence. 
Particles are described as material points which are dis- 
placed both by inertial forces (pressure-gradient force, 
added mass) and viscous forces (Stokes drag). Additional 
physical effects - as lift force, history force, buoyancy or 
finite-size and finite-Reynolds corrections - which may 
becomes important for the cases of light and/or large 
particles {[3 St > 1), are here neglected for simplicity. 

The particle dynamics then reads [18, 19] (see also [20]) 



dtx 



d t v = f3D t u + T l (u-v), (1) 



where x , v denote the particle position and velocity, re- 
spectively and dt = dt + v ■ V the time derivative along 
the particle path. The incompressible fluid velocity u 
evolves according to the Navier-Stokes equations 



D t w = d t u + u ■ Vu = —Wp/pj + isAu + f . 



(2) 



where p denotes the pressure and / an external forc- 
ing injecting energy at a rate e = (u ■ f). Eq. (2) is 
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evolved by means of a 2/3-dealiased pseudospectral code 
with a second order Adams-Bashforth time integrator. 
The fluid velocity at particle position is evaluated by 
means of a three-linear interpolation. Simulations have 
been performed in a cubic box of side L = 2?t with peri- 
odic boundary conditions, and by using N 3 = 128 3 and 
512 3 mesh points (reaching Taylor Reynolds numbers 
Rca = 75 and 180). The respective Taylor length scales 
A = y/ (u 2 ,) / (d x u^.) are A = 13r] and 21-q. The parameter 
space /3 x StE [0:3] x [0:4] is sampled with 504 (j3,St)- 
points with TV = 10 5 particles per type in the former case 
and 64 (optimally chosen by means of a MonteCarlo al- 
location scheme based on lower resolution results) with 
N = 1.6 • 10 6 in the latter. Given the small Re\ depen- 
dence, we will report here mostly results from Rc^ = 75 
as for that case we have a more complete sampling of the 
parameter space (j3,St). 

A requirement for any segregation indicator is to re- 
sult in zero segregation length for any two statistically 
independent distributions of particles coming from the 
same class of particles, at least in the limit of infinitely 
many particles. If the observation scale is too small and 
the number of particle finite, even independent particle 
realizations of the same class of particles will artificially 
appear to be segregated. Therefore, the definition of seg- 
regation strictly requires to indicate the observation scale 
r, and it will be sensitive to the particle number. 

However, we aim at a robust observable. Classical and 
natural observables, as e.g. the minimal distance between 
different type of particles strongly depends on particle 
number and hence are not robust. Harmonic averages 
of particle distances could be sensitive to small scales 
and not be spoiled by the large scales, but the choice of 
the weight exponent is rather arbitrary. The use of a 
density correlation function (p\(x)p2{x + r)) [17], though 
possible, requires to introduce a coarse-graining scale (to 
define the densities) which may quantitatively affect the 
estimate. The mixed pair correlation function, or mixed 
radial distribution function [23] , is not bounded at small- 
scales for clustered distributions of point-particles. 

Our approach is inspired by Kolmogorov's distance 
measure between two distributions [21] and is based on 
particle densities coarse-grained over a scale r, which can 
be understood as resolution of a magnifying glass used to 
look at the segregation. The whole volume L 3 is parti- 
tioned into Ai(r) = (L/r) 3 cubes. We then define the 
following segregation indicator: 

M(r) 

S aua2 (r) = N J2 K 1 -<1 ■ ( 3 ) 

The subscripts ot\ and a.2 index the particle parameters, 
i.e., a\ = (f3i,Sti) and ct2 = (/?2, ££2), N a is the total 
number of particles of a-type, while nf that of particles 
contained in each cube i. The case ol\ = 012 should be 
considered as taking independent realizations of the par- 




FIG. 2: S ai , a2 (r) vs. r for ai = (J3i = 0, Sti = 1.1) (heavy 
type) and ai — {P2 = 3,5*42 = 1.1) (light type) for differ- 
ent particle numbers N (= N ai = N a2 ): (o) N = 10 5 , (x) 
N = 10 4 and (+) N = 10 3 . Dashed and dotted lines refer 
to homogeneously distributed particles samples at various N. 
The expected Poisson scaling behavior, S ai ,a 2 oc r~ 3 ' 2 , is 
also reported. Inset: r seg , defined by S ai , a2 (r seg ) — 1/2, as 
a function of N for both heavy vs. light particles case and 
Poissonian samples. While for the former r 3eg saturates as N 
increases, for the latter r seg goes to zero like r seg ,h <x iV -1 '' 3 , 
as expected for Possonian samples (see text for details). 

ticle distribution (otherwise S ai ,ai =0 trivially) so that 
S a ^ a sets the minimum detectable segregation degree. 

Let us first discuss the limiting cases of S ai a2 (r). 
First, it can vary in the range [0,1]. S ail0l2 (r) = 1 
means that the two distributions arc not overlapping 
when looked at resolution r. For small enough scales, 
i.e., r -C l/p 1 / 3 (which is the mean distance of two par- 
ticles with p = N/L 3 the particle number density) this 
holds for any realization, therefore lim r ^o S ai ^ a2 = 1. 
On the contrary lim,.^^ S ai _ ct2 = as the total number 
of particles of the two species is globally identical (as as- 
sumed here). These limiting cases are observed in Fig. 2. 
Clearly, S aiy0l2 is a meaningful indicator of segregation 
only if it does not depend too severely on the particle 
number N . Indeed, in Figure 2, S ai . a . 2 (computed for the 
red and blue distribution of the central panel of Fig. 1) 
shows only a very weak TV-dependence at sufficiently 
large N. This is in contrast with the behavior of (3) for 
two independent and homogeneously distributed particle 
realizations (also shown in Fig. 2). The latter case can be 
easily understood recognizing that in each box of side r, 
m is a Poisson random variable, so that we can estimate 
rii w pr 3 ± yj pr 3 , where the two terms come from the av- 
erage and the fluctuation contributions, respectively. In 
eq. (3) the average cancels and, summing the fluctuations 
over all the (L/r) 3 cells, one has the order of magnitude 
estimate S(r) ~ [L/r) 3 N 1 / 2 ^ / L) 3 / 2 = N x ' 2 (r j L)' 3 ' 2 , 
explaining both the observed scaling behavior and the 
strong dependence on N. 

The segregation indicator allows us to extract the de- 
sired segregation length scale r seg . This can be done 
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FIG. 3: S ai , a2 (r) for ai = (J3i = 0, Sir = 1.1) and a 2 = 
{(52, St 2 = Sti), i.e., a heavy particle with /3i = and a given 
St vs. those having the same St but different densities 02- 
From bottom to top: fa = 0,0.5, 1, 1.5,3. Inset: r 3eg vs. fa, 
with r seg defined as in Fig 2, i.e., S aiiC , 2 (r seg ) = 1/2. The 
straight dashed line shows r seB) h ~ 4r/ ~ 0.3A. 

by fixing an arbitrary threshold value for S; we em- 
ployed S ai ^a 2 (r seg ) = 1/2 (see Fig. 2). With this def- 
inition, as shown in the inset of Fig. 2, for truly seg- 
regated (heavy vs. light) samples r seg saturates with 
increasing N. This does not hold for uniformly dis- 
tributed (non-segregated) particles. In the latter case 
r seg essentially coincides with the interparticle distance 
r seg = r seg ,h « l/p 1 / 3 = L/N 1 / 3 , as also seen from the 
inset of Fig 2. The behavior of r seg encompasses the fact 
that for a finite number of particles N a natural cut- 
off distance exists (the mean inter-particle distance) al- 
though we know theoretically that limTv^oo r seg (N) = 0. 
Hence r se g,h(N) can be interpreted as the accuracy in 
estimating r seg given a finite particle number N . 

We now calculate r seg for a pair of two different par- 
ticle classes to quantify their mutual segregation. Figure 
3 displays S aia2 (r) for distributions composed of heavy 
particles with Pi = and Sti = 0(1) and particles with 
the same St(= St 2 — Sti) but different densities. As 
one can see, segregation increases with the density dif- 
ference, but r seg w r se(Jt h for P 2 < 0.5, meaning that 
heavy enough particles basically all visit the same lo- 
cations in the flow, irrespective of their exact density: 
They tend to avoid vortical regions [8]. A sensitive in- 
crease of r seg is observed for (3 2 > 1 (Fig. 3 inset) and as 
expected the maximal segregation length is obtained for 
bubbles, i.e., particles with density ratio (3 = 3, where 
r seg 25?7 ~ 1.9A. For the same case, St = 1.1, (3 = 
vs. /3 = 3, at Re x = 180 we find r seg w 29ry « 1.4A. 

Thanks to the large number of particle types in our 
database we can extend the study of the segregation 
length to a wide range of physical parameters. In Fig. 4 
we show the value of the segregation length by fixing 
ai = ((3i, Sti) = (0, 1.1) (left), the red particles of central 
panel of Fig. 1 or the blue ones by fixing ai = (3, 1.1) 
(right) and varying a 2 = {P 2 , ^2) for the second kind of 




FIG. 4: r 3 e g between particle distributions with (3 = 0, St = 1.1 
(left) and /3 = 3, St — 1.1 (right) vs. distributions with generic 
(3, St. (•) indicates the reference particle type, and (x) the 
location of the maximal segregation length r a ™ g ax \ The solid 
contour line, traced at r seg = r aeg j- l = L/N 1 ^, sets the sen- 
sitivity level to distinguish between segregated and unseg- 
regated particle distributions. Dashed and dotted lines are 
drawn at r seg = n ■ r aeg> h, with n = 2, . . . , 6. The color scale 
codes the value of r aeg in units of the Kolmogorov length, r\. 

particles. The emerging picture is as follows. Particle 
class pairs with Sti w St 2 and (3i w (3 2 have a segregation 
length close to the interparticle distance r seg ^ and are 
unsegregatcd, while as soon as the Stokes number or the 
density difference become larger, r seg >r seg j l . The max- 
imal segregation length (r^g^ «27?7«2.1A) is roughly 
twice the Taylor microscale and is realized for particles 
with large density difference (3\ = and (3 2 >1 (or (3i = 3 
and (3 2 <1). These results thus confirm those of Fig. 3. It 
is interesting to note that heavy couples, (3i, (3 2 <1, segre- 
gate less than light ones, (3i,(3 2 >l, which arc thus much 
more sensitive to small variations of density and/or re- 
sponse times. The correlation between position and flow 
structure is thus much stronger for light particles. 

We next systematically study what happens when fix- 
ing St or (3. In the first four panels (from left) of Fig. 5, 
r aeg is shown as a function of (3i and (3 2 , where we fixed 
the Stokes numbers, Sti = St 2 = St= (0.31, 0.6, 1.1, 4.1). 
Close to the diagonal (3i ps (3 2 it is r seg < r seg _h, i.e., 
similar particles are poorly segregated. Outside these 
regions r seg > r se g,h an d segregation is above the ac- 
curacy threshold r se g,fe- Two observations arc in order. 
First, the strongest segregation is present for the case 
of St ~ 1 meaning that response times of the order of 
the Kolmogorov time are best suited to generate strong 
correlations between flow structures and particle posi- 
tions and consequently segregation. Further, segregation 
is stronger for couples composed by very heavy ( (3 ~ 0) 
and very light ((3 ~ 3) particles. Therefore, light and 
heavy particles with St ~ 1 display the strongest cluster- 
ization. Second, the numerical value of the segregation 
length saturates to a constant value ~ 2A, strongly in- 
dicating that we are measuring an intrinsic property of 
the underlying turbulent flow emerging when particles 
are strongly clusterized. 

The two rightmost panels of Fig. 5 display cuts done 
by fixing the density, namely (3i = (3 2 — (3= 1.5 (rcsp. =3). 
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FIG. 5: From left to right, segregation length r 3eg between distribution of particles with St = (0.31,0.6,1.1,4.1) vs. the 
densities /3i,/?2 and (last two panels) for particle class pairs having the same densities /3i — (3% — 1-5 (resp. 3) and different St. 
The color scale and the contour lines are, as in Fig. 4, at r seg =n ■ r seg> h, with n = 2, . . . , 4. 



For (3 < 1.5, i.e. relatively heavy or weakly light parti- 
cles, there is only a slight tendency toward segregation at 
varying the Stokes number. This means that even if the 
particles form clusters, such clusters are not too sensitive 
against variation of St. For very light particles, (3 = 3, the 
situation is different when comparing the case of St ~ 1 
with a St ~ case. As expected for 5f~0, though parti- 
cles are light, they distribute almost uniformly (they are 
not too far from the tracer limit) while for St>l they 
are strongly correlated with the vortex filaments which 
are unevenly distributed in the flow [22] . 

In conclusion, we introduced an indicator able to quan- 
tify the segregation degree and allowing to define a seg- 
regation length scale r seg between different classes of 
particles, which follow simplified dynamical equations in 
isotropic turbulence. The extracted information is in line 
with the intuitive idea of expulsion/entrapment of par- 
ticles due to vortical structures which is now on a more 
quantitative ground. The maximal segregation length, 
for instance for heavy particles ((3 = 0, St =1.1), is ob- 
tained with bubbles with slightly larger Stokes (54 ~ 1.4) 
(Fig. 4 left); it measure r^eg 3 ^ ~27?? ~2.1A at Re\ = 75. 
At Re x = 180 we get similarly: ri™ x) - 48?y ~ 2.3A. 
Therefore, rj^ g is about twice the Taylor length. 

Important areas of application of the developed meth- 
ods go far beyond homogeneous isotropic turbulence in- 
cluding, e.g., heterogeneous catalysis or flotation, where 
one is interested in the collision probability of argon bub- 
bles and solid contaminations in turbulent liquid steel 
[16]. In this context our finding that r^g ~ 2A sug- 
gests that the cleaning could become less efficient at high 
Rex, as bubbles and particles then become more seg- 
regated. Hocver, quantitative statement requires bet- 
ter models for both the flow geometry and the effective 
particle force. Another example for the application of 
the suggested method is the formation of rain drops at 
solid nuclei in clouds [17], a mechanism which is crucial 
to develop models for rain initiation [6]. Disregarding 
particle segregation in all these examples would lead to 
estimates of the collision rate which could be orders of 
magnitude off. We have already observed that the point- 
particle model adopted in our study neglects some hydro- 
dynamical effects which corresponds to additional forces 



in the particle equation. Such forces, when included with 
proper modeling, might smooth the intensity of segre- 
gation which, as we have shown, is mainly due to the 
relative strength of inertial forces. 

Finally, we stress that the introduced segregation in- 
dicator can be employed in all phenomena involving dif- 
ferent classes of segregating objects. Provided one knows 
the position of all the objects, no prior knowledge on the 
physical mechanism which leads to segregation is needed. 

We thank J. Bee and L. Biferale for fruitful discussions 
and R. Pasmanter for prompting Kolmogorov's distance 
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